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We calculate the back-reaction of long wavelength cosmological perturbations on a general relativistic 
measure of the local expansion rate of the Universe. Specifically, we consider a cosmological model 
in which matter is described by two scalar matter fields, one being the inflaton and the other one 
representing a matter field which is used as a clock. We analyze back-reaction in a phase of inflaton- 
driven slow-roll inflation, and find that the leading infrared back-reaction terms contributing to the 
evolution of the expansion rate do not vanish when measured at a fixed value of the clock field. We 
also analyze the back-reaction of entropy modes in a specific cosmological model with negative square 
mass for the entropy field and find that back-reaction can become significant. Our work provides 
evidence that, in general, the back-reaction of infrared fluctuations could be locally observable. 



I. INTRODUCTION 

Because of the nonlinear nature of the Einstein equa- 
tions, the ansatz for metric and matter used all the time 
in early Universe cosmology, namely a homogeneous and 
isotropic background space-time plus fluctuations which 
average to zero on the spatial hypersurfaces (set by the 
background cosmology) and which obey the linear fluctu- 
ation equations, does not obey the Einstein equations at 
second order in the expansion parameter (the fractional 
amplitude of the linear fluctuations). Terms which have 
to be added to the metric such that the Einstein equa- 
tions are satisfied to second order are called back-reaction 
terms. 

The back-reaction of short wavelength gravitational 
waves on an expanding Friedmann-Robertson- Walker 
(FRW) cosmology is a well-established subject [1]. We 
insert the ansatz for the metric consisting of linear fluc- 
tuations Sg^y about the FRW background 



9fiu = 9% + Sg^ 
into the Einstein equations 



Guv — T, 



(1) 



(2) 



G M „ and denoting the Einstein tensor and the energy- 
momentum tensor, respectively, expand to second order, 
use the fact that the linear fluctuation equations are sat- 
isfied to cancel the linear terms, and move all quadratic 
terms to the right hand side of the equation. These 
quadratic terms are interpreted as an effective energy- 
momentum tensor t^ v ( "pseudotensor" ) for fluctuations. 
The resulting equation 



G^(g br )=T^ br ) + T f . 



br \ 



(3) 



(up to terms of third and higher order) can be inter- 
preted in the sense that the fluctuations carry energy 



and pressure and hence effect the background space-time 
in which they are defined. Note, that in order for this 
equation to be satisfied, we need to add quadratic cor- 
rection terms ( "back-reaction terms" ) to the metric and 
the matter (which explains the superscripts in the above 
equation). In particular, the homogeneous component of 
the metric obtains a correction term of quadratic order 
which can be extracted from (3) by taking spatial aver- 
ages (see [2,3] for more details). 

Since, in the context of inflationary and post- 
inflationary cosmology, the scalar metric fluctuations 
(fluctuations coupled to energy density and pressure per- 
turbations) are believed to dominate over the effects of 
gravitational waves, it is of great interest to study the 
back-reaction of these cosmological perturbations. Fur- 
thermore, in inflationary cosmology the phase space of in- 
frared modes (defined as modes with wavelength greater 
than the Hubble radius) is growing over time, whereas 
the phase space of ultraviolet modes does not grow, and 
since the amplitude of the associated metric fluctuations 
of these infrared modes does not decrease in time (see 
e.g. [4] for a comprehensive overview of the theory of cos- 
mological fluctuations and [5] for a recent introductory 
overview) , the back-reaction of these infrared modes may 
be very important. 

Rather recently, the effective energy-momentum ten- 
sor formalism of [1] was generalized to study the back- 
reaction of infrared modes on the spatially homogeneous 
component of the metric [6,2]. The surprising result 
which was found is that in a slow-roll inflationary back- 
ground, the infrared contribution to t^ v has the form 
of a negative cosmological constant, with the absolute 
value of the energy density increasing in time as the phase 
space of infrared modes grows and the amplitude of the 
gravitational potential remains large (in fact the ampli- 
tude is slowly increasing). A similar result was found in 
an earlier study by Tsamis and Woodard [7], in which 
the back-reaction of long wavelength gravitational waves 
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in a de Sitter background was studied *. The results 
of [6,2] were confirmed in [21] using very different tech- 
niques. On the basis of these calculations, it was argued 
[6,2,21] that gravitational back-reaction in scalar field- 
driven inflationary models, calculated up to quadratic 
order in perturbations and to leading order in the long 
wavelength expansion and in the slow roll approximation 
could decrease the expansion rate of the universe and po- 
tentially solve the cosmological constant problem [22,23] 
(see also [7,24] for similar suggestions in the context of 
back-reaction studies of gravitational waves, and [25] for 
earlier ideas concerning the instability of dc Sitter space) . 

However, as was emphasized by Unruh [26] \ the ap- 
proach of [6,2,7] has several shortcomings. First of all, 
due to the nonlinear nature of the Einstein equations, 
calculating an "observable" from the spatially averaged 
metric will not in general give the same result as calculat- 
ing the spatially averaged value of the observable. More 
importantly, the spatially averaged metric is not a local 
physical observable. An improved analysis of gravita- 
tional back-reaction on the locally measurable expansion 
rate of the Universe starts with identifying a local phys- 
ical variable which describes this expansion rate, then 
calculating the back-reaction of cosmological perturba- 
tions on this quantity, and only at the end taking an 
expectation value of the result. It is important to fix the 
hypersurface of averaging by a clear physical prescription 
in order to remove the possibility of being misled by co- 
ordinate artifacts which are not physically measurable. 
Such an approach was first suggested by Abramo and 
Woodard in [27]. 

In a first paper [28], we considered a simple variable 
describing the local expansion rate and calculated this 
variable to second order in the metric fluctuations in a 
model with a single matter field, the inflaton field if. 
We then evaluated the leading contributions of infrared 
fluctuation modes to this variable. When evaluated at a 
fixed background time t, we found that the back-reaction 
effects were present. However, when evaluated at a fixed 
value of the matter field if, the only physical clock avail- 
able in this simple system, the dependence of the expan- 
sion rate on the clock time took on exactly the same form 
as in an unperturbed background. Thus, the leading in- 
frared back-reaction terms had no locally measurable ef- 
fect in this system. Similar conclusions were reached in 
[29]. 



However, the fact that the infrared back-reaction terms 
did not vanish when evaluated at a fixed background time 
leads us to the conjecture that they will have a locally 
measurable effect provided we have a clock different from 
the inflaton field ip itself The simplest way to introduce 
an independent physical clock is to add a second scalar 
matter field \ to the system. This field \ represents our 
observed matter fields. 

In this paper, we evaluate our variable describing the 
local expansion rate of the Universe to second order in 
the fluctuations in the slow-rolling phase of an inflation- 
ary Universe with two matter field, the inflaton ip which 
dominates the energy-momentum tensor and a spectator 
(clock) field x- We find that, in general, the leading in- 
frared terms do not vanish when evaluated at a fixed 
value of the clock field \. Thus, in this system infrared 
back-reaction is for real We also discuss some subteties 
concerning this conclusion. The physical effects of this 
infrared back-reaction remain to be investigated [31]. 

The outline of this article is as follows: 

In the following section we review the construction 
of the variable 9 which we propose to use to describe 
the local cosmological expansion rate, and we summarize 
the results of our previous paper [28] in which this vari- 
able was computed in a cosmological model with a single 
scalar matter field. 

In Section 3, we derive the expression for 6 in a model 
with two matter fields in terms of the values of the met- 
ric fluctuations (more specifically, in terms of the gravi- 
tational potential $ which characterizes the fluctuations 
in the longitudinal gauge). 

In Section 4, we derive an expression for the change of 
G over a short time during which the clock field \ changes 
by a small amount, in terms of physically measurable 
quantities, we investigate the behavior of back-reaction 
in scenario where adiabatic modes are dominant and \ 
is a rapidly oscillating scalar filed where the local energy 
density of x lSi our clock. We show that when evaluated 
during the slow-roll phase of an inflationary Universe, 
there are non- vanishing back-reaction contributions from 
infrared modes. 

Then, in section 5, we consider a two matter field 
model motivated by hybrid inflation, where it is expected 
that the entropy fluctuations are large. We calculate the 
back-reaction effect in terms of proper time, as a clock 
and we show that back-reaction of entropy modes could 
be significant. 



"The back-reaction of small-scale (i.e. smaller than the Hub- 
ble radius) cosmological perturbations has been considered in 
[8-15], and in [16,17] in the context of Newtonian cosmo- 
logical perturbation theory. Note also that Nambu [18-20] 
is developing a program to compute back-reaction effects on 
the spatially averaged metric using the renormalization group 
method. 

T We are also grateful to Alan Guth, Andrei Linde and Alex 
Vilenkin for detailed private discussions on these points. 



*To explain the effects of back-reaction in heuristic lan- 
guage one needs to introduce a second matter field: the back- 
reaction effect is the fact that the expansion rate is different 
in the presence of fluctuations compared to the result in the 
absence of fluctuations when measured at the same value of 
the temperature of the CMB. 

§ Another system in which infrared back-reaction is for real 
was discussed in [30]. 
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In the final section we summarize and discuss our re- 
sults and point out issues which remain to be resolved. 

II. LOCAL BACK-REACTION IN SINGLE FIELD 
MODELS 

We begin by reviewing the construction of the variable 
measuring the local cosmological expansion rate: We 
consider the velocity four-vector field u a tangential to a 
family of world lines of a matter fluid in a general space- 
time. This four-vector is normalized such that 

u a u a = 1 , (4) 

where a runs over the space-time indices. In terms of 
this four-vector, the local expansion rate 6 is defined by 

Q = u« a . (5) 

It represents the local expansion rate of the tangential 
surfaces orthogonal to the fluid flow. 

For a homogeneous Universe with scale factor a(t), the 
Hubble expansion rate H is related to 6 via 

H = d t a/a=^@. (6) 

Thus, for a Universe with fluctuations, we use dtS/S = 
|0 to define the "locally measured" Hubble expansion 
rate. This variable is a mathematically simpler object 
than the one recently introduced in [27], which involves 
the integral along the past light cone of the observation 
point of a function which itself is rather difficult to cal- 
culate. 

In [28], we considered a theory with a single matter 
field <p and evaluated 6 up to second order in the met- 
ric fluctuation variable $ (defined later in (16)). When 
expressed as a function of physical time t, and dropping 
all spatial gradient terms, the result obtained is 

Q 

9 = 3-(l - $ + -$ 2 ) - 3($ + $$) (7) 
a £ 

Thus, at first glance it looks like the local expansion rate 
is different in the presence (6 ^ 0) than in the absence 
(6 = 0) of fluctuations. However, the above conclusion is 
premature since the background time t is not a physical 
quantity. 

In order to draw conclusions about the local measura- 
bility of the back-reaction of infrared fluctuation modes, 
it is necessary to express the final result in terms of a lo- 
cally measurable quantity, a clock. The only such clock 
available in the single matter field model considered in 
[28] is the inflaton field itself. Thus, we must express 
the metric fluctuations $ in terms of the matter fields. 
This can be done by making use of the Einstein con- 
straint equations. Neglecting spatial derivatives in these 
equations, one can show (see also [32]) that, also setting 



4> = (which is approximately the case for each long 
wavelength fluctuation mode in a slow-roll inflationary 
model [4]), the final result for the local expansion rate 
becomes 

e = Vs^v^p) (8) 

which as a function of ip is the same as the relation for 
an unperturbed background. Thus, in the single matter 
field model the leading back-reaction terms from infrared 
modes are not locally measurable, in agreement with the 
conjectures on [26] ** 

In the presence of a second matter field \-> it is possible 
to use this field as a clock. In fact, using this field to 
model the temperature of the CMB, it is physically more 
sensible to use \ as the clock rather than the inflaton 
field itself. The fact that in the single field matter model 
the dominant infrared back-reaction terms do not vanish 
when evaluated at a fixed background time t leads us to 
conjecture that they will not vanish when evaluated at 
a fixed value of the new clock field \ in the two matter 
field model. 

If we are able to verify the above conjecture, we 
would have established a close analogy (already suggested 
in [28]) with the analysis of the parametric amplifica- 
tion of super-Hubble-scale cosmological fluctuations dur- 
ing inflationary reheating. From the point of view of 
the background space-time coordinates, it appears [34] 
that the parametric amplification of matter fluctuations 
on super-Hubble scales in an unperturbed cosmological 
background (see e.g. [35,36] for a discussion of parametric 
resonance during reheating) would imply the parametric 
amplification of the cosmological fluctuations on these 
scales. However, it can be shown that in single field 
models physical observables measuring the amplitude 
of cosmological fluctuations do not feel any resonance 
[37-39,32]. In contrast, in two field models of inflation 
there is [40,41] parametric amplification of super-Hubblc- 
scale cosmological fluctuations. In this case, there is a 
fluctuation mode corresponding to entropy fluctuations 
which cannot locally be gauged away. This mode is (in 
certain theories) parametrically amplified during reheat- 
ing, and in turn drives the parametric resonance of the 
super-Hubble scale curvature fluctuations. 

III. DERIVING THE EXPANSION RATE FOR 
MODELS WITH TWO MATTER FIELDS 

The procedure for calculating will be as follows: 
First, we determine the velocity four- vector field u a for 



"Note, however, that there still could be back-reaction of 
infrared terms which can change the coarse of evolution of ip, 
due to the fact that new modes are continuously exiting the 
Hubble radius. These effects are calculated in [33]. 
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the two field model. Then, we use the Einstein equations 
to express u a in terms of the metric perturbations. Tak- 
ing the relative amplitude of the metric fluctuations as 
the expansion parameter, we then calculate 6, our local 
measure of the Hubble expansion rate, to second order. 
After evaluating the result on a physically determined 
hypersurface we can then study the back-reaction of cos- 
mological fluctuations on the locally measured Hubble 
expansion rate. We will focus on the leading infrared 
contributions to back-reaction, the terms found to dom- 
inate the back-reaction effects in [2,6,21]. 

We consider a model with two scalar matter fields, the 
inflaton ip and the clock field \'- 



(9) 
(10) 



We describe the inflaton as a perfect fluid with energy 
density p and pressure P, considering 



(p + P)u^u v - Pg^ v = d^pdvp - C .g^ v 

which can be justified by taking 

P = C and 
1 



P = 2 d ^ d ^ + X) , 



which leads to 



itp = Ad^p 
A = (d v pd v p) 



-1/2 



(11) 

(12) 
(13) 

(14) 
(15) 



Starting from the expression for the metric to linear 
order in the fluctuations <3> (see [4] for a detailed review), 
we determine the velocity four-vector field u a to second 
order, the order required to analyze the leading infrared 
terms in the back-reaction to quadratic order. In linear 
perturbation theory, and in the case of simple forms of 
matter for which there is to linear order no anisotropic 
stress, the metric (in longitudinal gauge) can be written 
as 

ds 2 = a(r/) 2 ((l + 2*)rfr; 2 - (1 - 2$)j ij dx i dx j ) , (16) 



lij =S ij [l + -lC(x 2 +y 2 + z 2 )} 



(17) 



where JC = 0,1,-1 depending on whether the three- 
dimensional space corresponding to the hypersurface t = 
const, is flat, closed or open. In this paper we will take 
it to be zero in order to simplify the calculations. The 
time variable rj appearing in (16) is conformal time and is 
related to the coordinate time t via di] = a~ 1 dt. For the 
forms of matter considered here, * = $ at linear order 



tt . The field $ is then called the relativistic gravitational 
potential. 

In order to obtain the complete result for gravitational 
back-reaction, we consider the Einstein equations for a 
perfect fluid and additional scalar field \. 

G>„ = (P + p)u^u v - Pg^ + d^xduX - £'dnv i 18 ) 

(in units in which 8irG = 1), which, since = —R, will 
yield 



-R = p-3P-d» X d^X- 



(19) 



We can define a new tensor called G'^ v and a new scalar 



Rf by, 



G' 



R' = R-d» X duX 



(20) 
(21) 



Then, the equations become similar to the case of a single 
scalar field (see [28]), and p satisfies the equation, 



p = vTG'^u* 



(22) 



which leads to an equation that can be solved perturba- 
tively for Ui and w in parallel with equation (4) to any 
desired order: 



Gf - -u"G^«"« tij + R'u° Ui . 



(23) 



For an unperturbed Robertson- Walker metric, the 
four- velocity field u in comoving coordinates would be 



u" = (1,0,0,0), 



(24) 



which can be substituted in cquation(23) as zeroth order 
solution in an expansion in powers of <!> to calculate u 1 
up to second order, 



, ( 2) G M - d° X d i X 



(25) 



where G 0t is proportional to d l Q. Using equation (4) we 
can derive the expression for the time component of u a 
in terms of $: 

u°(r]) = a - 1 (l-$ + ^$ 2 ) + iauV. (26) 

Now that we have all components of u a , we need to 
take the covariant derivative of this vector and retain all 
$ dependence up to second order. During inflation, fluc- 
tuations which are generated on sub-Hubble scales early 
on during the inflationary phase are red-shifted to scales 



^Even if we did not make the assumption VP = 3>, it turns 
out that at second order all infrared terms depending oii*-$ 
will drop out. 
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much larger than the Hubble radius. Thus, in this con- 
text it is of great interest to consider the back-reaction 
of infrared modes. To investigate this effect we neglect 
any spatial derivative term contributing to the local ex- 
pansion rate 6. Since 



V9 



(27) 



and, according to cquation(25), spatial components of 
u are proportional to spatial derivatives of $ and x, we 
conclude that, neglecting the spatial derivative terms, the 
expression for 6 reduces to 



6 



1 



A straightforward calculation yields the result 
9 = 3-(l - $ + ^3> 2 ) - 3(4 + $$) . 



(28) 



(29) 



Where a dot denotes differentiation with respect to the 
coordinate time t. Since H = a/ a, we can immediately 
read off the extra terms contributing to the local expan- 
sion rate which result from the presence of cosmological 
fluctuations. Hence, it follows that if evaluated at a fixed 
coordinate time, infrared modes on average lead to cor- 
rections in the expansion rate compared to what would 
be obtained at the same coordinate time in the absence of 
metric fluctuations. Whether this is a physically measur- 
able effect will be discussed in more depth in the following 
section. 



IV. EXPANSION RATE IN TERMS OF A 
SCALAR FIELD AS OBSERVABLE 

The expression (29) gives 9 in terms of the metric 
fluctuations. However, to determine whether this local 
effect of gravitational back-reaction of cosmological fluc- 
tuations is physically measurable, we need to express 6 
in terms of a physical clock field, not just in terms of 
the background time. As was mentioned before, we will 
use x as a clock. Hence, to obtain results with physical 
meaning, we must evaluate 6 on a surface of constant \ 
and not constant t. 

By dropping the gradient terms from the G 00 Einstein 
equations, it can be shown that 



H 



VI + 2$ 



+ 2t ] 
yoo £ yoo 



(30) 



where (1 + 2$) is simply g ao , the time-time component 
of metric. If we expand the left hand side in terms of <f> 
we get the result, 

(31) 



Substituting this result into equation (29) and neglecting 
terms containing $ (which are sub-dominant compared 
to other terms in the infrared limit as long as the equation 
of state of the background is not changing significantly 
over time) leads to Friedmann-like equation for the local 
expansion rate 



* yoo yoo 



(32) 



At this point, it appears at first glance that the expansion 
rate for a fixed value of the clock field manifestly depends 
on the value of the gravitational potential $ which enters 
via g 00 , and that thus the leading infrared back- reaction 
effects are locally measurable. However, we are not yet 
ready to reach this conclusion since the right hand side 
of (32) still contains derivatives with respect to the back- 
ground coordinate time which should be re-expressed in 
terms of local physical observables. 

To describe the behavior of 9 on surfaces of constant 
X, we have to express the right hand side of (32) in terms 
of x arj d of the initial values tfi n and Xin of the matter 
fields. This involves solving the Klein-Gordon equations 
for both scalar fields: 



H 



1 



V, v (<p, X ) + 3 — <f + — </3 = 

Qoo Qoo 

v x {ip, x ) + i— x + — x = o. 



(33) 
(34) 



Since finding exact and analytical solutions of this cou- 
pled system of equations is not possible for a general 
form of the potential, we solved these equations for an 
infinitesimal period of time St. In other words, we used 
solutions linear in time to investigate the dynamic behav- 
ior of 9. We start with initial values of ip in (x), <j>i n , Xin 
at ti n on a surface with a constant value of Xin, let the 
system evolve, and calculate the final value of these pa- 
rameters on a surface with constant value of Xf- The 
infinitesimal evolution equations are 



St = Xf . Xm 

Xin 

<Pf = <fiinSt + (pi n 
Iff = (finSr) + ip in 
Xf = XtnSt + Xtn ■ 



(35) 

(36) 
(37) 
(38) 



Notice that we do not assign a uniform time to all points 
on the initial hypersurface, and that, for now, we just 
solve equations locally. Inserting the results from the 
above equations into the general expression (32) for 9/ 
leads to 



e / = e ir , 



3(X/ ~ Xin) 
^Xin \/Qoo 



($n + Xin) . 



(39) 



where the square root of g 00 is evaluated on the initial hy- 
persurface. We are still not finished, because the second 
term on the right hand side still contains derivatives with 
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respect to the coordinate time. So we rewrite the equa- 
tion (39) using the proper time instead of the coordinate 
time (<9t = Jgoodt), for gauge invariance purposes, 



v(x) = 



Am 



(40) 



where 5\ — Xf ~ Xin- Equation (39) can be further 
simplified using Equations (30) and (32), and this enables 
us to express 6/ as a function of 0i„, <pi n , x and x'i 



Of = ®in- 



3 2 
2 { 3 



- V((p in ,Xin)) 



5x 



(41) 



This is the completely local result for the change in O 
which we were aiming to obtain. One may think that 
(41) implies there is no local back-reaction **. However, 
this conclusion is not correct because in the presence of 
cosmological fluctuations the trajectories of the x an d f 
fields start to deviate from the classical solution that are 
obtained from Equations (33) and (34). Thus, the values 
of ip' in and x'in f° r the next time steps will fluctuate on 
the Xin = const hypersurface, whereas in the absence of 
such fluctuations, the initial values at each step follow 
the classical trajectory. In other words, compared to the 
classical solution, there is a difference in the evolution of 
due to the different paths for each point of an inhomo- 
geneous universe in field space. Furthermore, there are 
other terms coming from the modes that are exiting the 
Hubble radius and are not taken in to account in the in- 
frared limit of Equation (41). Interpreted correctly, our 
result (41) thus shows that one expects that in a two mat- 
ter field model, infrared cosmological fluctuations have a 
locally observable effect if cosmological perturbations are 
excited. 

To clarify this point we investigate the behavior of 
back-reaction terms in two scenarios, first where adia- 
batic modes are dominant and then, in the next section, 
in a scenario where entropy modes play a major role. 
In the former case we assume that the x ne ld does not 
significantly contribute to the energy density of the uni- 
verse and only plays the role of a clock. We can consider 
X to be a rapidly oscillating scalar filed (to mimic the 
cosmic microwave background radiation in our Universe) 
and take the local energy density of x as our clock. For 
that purpose we work with the following potential 



**The completely local form of Equation (41) is easy to un- 
derstand from the mathematical point of view: by neglecting 
all of the spatial derivative terms, not only are we treating 
our equations completely locally, but we are also eliminating 
all interactions between neighboring points in space. Thus, 
the final result for 9 at a specific point in space can depend 
on nothing more than on the initial values of the fields ipi n 
and Gi n at that same point. In other words, the dynamics 
of the fields at each point of space behaves like a Friedmann 
Universe with the appropriate initial conditions. 



(42) 



where, by taking p ^> 0, the rate of oscillation of x will 
be much faster than the expansion rate, which enables 
us to solve Equation (34) using the WKB approximation 
and to calculate the value of p x — x' 2 /2 + V(x), the 
energy density of x, 



Px = P x ,mcxp(- / Qdr) . 

Jr in 



(43) 



For the inflaton field, we also take a standard quadratic 
potential of chaotic inflation, 



V(<p) = -mV- 



(44) 



To take into account the effects of the ultraviolet modes 
that become infrared as the inflationary comoving Hub- 
ble radius shrinks, we modify Equation (33) in the slow- 
roll regime to obtain 



0(^' + mV = C Q 5/2 C(x,T) 



(45) 



where C a (a standing for adiabatic) is a numerical con- 
stant and £ is a Gaussian random function with the co- 
variance 

(£(xi,ti)£(x 2 ,t 2 )) = S(n - T 2 )smc(aH\x 1 - x 2 |) . (46) 

In order to derive Equation (45) rigorously, one needs 
to coarse-grain the non-linear Einstein equations on the 
scale of the Hubble radius §§. The stochastic term on the 
right hand side of Equation 45 is due to the adiabatic 
fluctuation modes that are exiting the Hubble radius (see 
[33] for a detailed derivation). 

Now, using the fact that Equation (32) simplifies to 
= y/3V{<p) in the slow-roll approximation, one can 
solve Equation (45) perturbatively in £, (see [33] for de- 
tails) and calculate the expectation value of 0, 

~ (0) = m 2 (-r) + ^L(0?„ - m 6 (-r) 3 ), (47) 

where we have taken r to be negative during and to go 
to zero at the end of the slow-roll period. This result 
shows that back-reaction increases the expansion rate 
when proper time is used as the clock. As we shall see, 
this result changes when we use p x as a clock. To sim- 
plify the expression for the expansion rate when p x is 
used as a clock, we define a new parameter, (3 = In 

which is a measure of p x . Here, p fc is a constant and the 



^This is similar to the approach taken in deriving the equa- 
tions of the stochastic scenario of inflation [42]. A more rig- 
orous derivation which includes metric perturbations is given 
in [32]. 
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subscript fc refers to the time t. = 0, the end of the 
classic slow-roll inflation. Hence, we have 



-f 

Jo 



e(f)df. 



(48) 



Now, substituting 6 from Equation (47), in the above 
integral we can calculate r in terms of 0, 



vw 



to 



{1 



r 2 i Pi 3 

a ( m 4 3 2 + in 
2/3V6 3 M 3m 



)}■ (49) 



Substituting r back in Equation (47), we can drive in 
terms of 0: 



~ my/20 



1 %r 2 

^±m 3 (0f/ 2 



(50) 



Interestingly, the sign of the back-reaction term is neg- 
ative, whereas it was positive in Equation (47). Thus, 
it appears that when we measure local expansion rate 
using an oscillating field as our clock, unlike when we 
used proper time, back-reaction slows down the expan- 
sion rate. Describing the back-reaction in terms of is of 
more relevance to the inflationary predictions for today's 
universe, as is proportional to number of c-foldings of 
inflation which enters into the calculation of the power 
spectrum of metric fluctuations. 



V. BACK-REACTION OF ENTROPY MODES 

In this section we will consider a two matter field model 
motivated by hybrid inflation [43], where it is expected 
that the entropy fluctuations are large. In this section, 
due to the specific choice of the model that we inves- 
tigate, we calculate the back-reaction effect in terms of 
proper time, i.e. constant- t surfaces, as a clock instead 
of constant-x surfaces. This is an appropriate choice be- 
cause we only want to investigate the magnitude of this 
effect. In terms of this clock variable, Equation (41) takes 
the following form, 



Qf = 0,; 



3 2 
2 [ 3 



in i "Kin 

))5r 



(51) 



where we see that different trajectories in the (<p,x) 
space, which could belong to different points in an in- 
homogeneous universe, lead to different behaviors for 0. 
As we demonstrate below, this effect could be clearly seen 
when we consider the evolution of 6 in the long run. 

A simple model where we get growing entropy modes 
is hybrid inflation, where the secondary field, x> has a 
negative square mass. In the slow-roll regime, the clas- 
sical value of x vanishes, which is why we do not use it 
as our clock. Now, notice that the local evolution of 9 is 
given by 



Q f = e in -l f\x' 2 + v ,2 )dT. 



(52) 



If we can show that, due to the presence of entropy 
modes, x' 2 may dominate ip' 2 in Equation (52), we have 
demonstrated that back-reaction can dominate the evo- 
lution of 0. 

For the sake of clarity, let us investigate the behavior 
of x' 2 in a model where, 



V(<p,x) = ^mV - \^X 2 - 



(53) 



If we assume x is also slowly rolling, then, similar to what 
happens in Equation (45), the field equation (34) should 
be modified to take into account the entropy fluctuation 
modes, while we ignore such modifications (arising due 
to the adiabatic modes) in Equation (33). This leaves us 
with the following field equations 



C e 5 / 2 £(x,r), 



Qip' + m 2 tp = 0, 



(54) 
(55) 



where C e (e standing for entropy) is a numerical constant. 
Again, using ~ y/3V(<p), we can calculate x up to first 
order in £, 



C e m\-Tf 



-fj, 2 /m' 



/: 



(-T) 



~\3/2+ fJ , 2 /m' 



£(x, f)df, 



which yields 



4 + 2 A? T 



(56) 



(57) 



where /(r) is a correction due to the fluctuations of the 
Hubble size modes, and its amplitude decreases with 
time. For slow-roll conditions to be satisfied, we need 
to have r > — . Inserting this limit for the growing mode 
in (x' 2 ) and comparing it with homogeneous value of ip' 2 , 
which is the classical rate of change in 0, we get 



<x' 2 > 



.12 



477J 4 



4 o 
+ 2^ 2 TO 2 V 



4+2 M 2 /™ 2 



(58) 



at the end of slow- roll inflation. Since ifi n > 1, this rela- 
tion implies that if /i > to, then back-reaction dominates 
the inflation before its end (calculated in the absence of 
back-reaction), and the homogenous background solution 
loses its validity. 

This result can also be important for calculating the 
spectrum of fluctuation after inflation ends. We have 
also seen that since different points have different initial 
conditions, they will have different final expansion rates. 
To obtain a quantitative measure of the back-reaction 
of infrared modes, one needs to take an average of the 
expansion rate over the many initial Hubble patches that 
form our observable universe. 
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VI. CONCLUSIONS 

In this paper, we have studied the leading back- 
reaction effects of long wavelength cosmological fluctu- 
ations on a local observable 6 which measures the local 
expansion rate of the Universe. The observable gives the 
rate at which neighboring comoving observers separate 
and coincides with the usual definition of expansion in 
the context of the fluid approach to cosmology. We con- 
sidered models with two matter fields, an inflaton field ip 
which dominates the energy-momentum tensor and leads 
to slow-roll inflation, and a second matter field \ which 
was used either as a clock or it contributed to produc- 
tion of entropy modes. We investigated the behavior of 
back-reaction terms in two scenarios, first where adia- 
batic modes are dominant and then, in a scenario where 
entropy modes play a major role. In the former case to 
obtain locally measurable statements about the presence 
or absence of back-reaction effects, We considered \ to be 
a rapidly oscillating scalar filed (to mimic the cosmic mi- 
crowave background radiation in our Universe) and took 
the local energy density of \ as our clock. We evalu- 
ated and found that not only the back-reacton term 
does not vanish but the sign of the back-reaction term is 
negative as well. In the second scenario, we analyzed the 
back-reaction of entropy modes in a specific cosmologi- 
cal model with negative square mass for the entropy field 
and found that back-reaction can become significant and 
dominate the evolution of O. 

Our main result is that the leading infrared terms, the 
terms which dominate the effects discussed in [2] and 
[21], arc non- vanishing and physically measurable if we 
use the second scalar field \ as a clock. The effect can be 
particularly large if the entropy mode of the cosmological 
fluctuation field is excited. This confirms the conjectures 
made in our previous work [28] and is in good agreement 
with the physical intuition gained from the study of the 
parametric amplification of supcr-Hubble length cosmo- 
logical fluctuations during reheating [40,41]. 

In this paper, we have thus shown that infrared back- 
reaction is "for real" in two matter field models. How- 
ever, how important the back-reaction is remains to be 
investigated. Results will be reported in a subsequent 
publication [31]. 
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